
CYCLOTRON HARMONIC WAVE PROPAGATION AND 
INSTABILITIES: II. OBLIQUE PROPAGATION 

by 

J. A. Tataronis and F. W. Crawford 



AEC Contract AT (04- 3) 32 6 
(Project Agreement No. 1) 

and 

NASA Research Grant NGR 05-020-077 


SU-IPR Report No. 326 
June 1969 


Institute for Plasma Research 
Stanford University 
Stanford, California 



CYCLOTRON HARMONIC WAVE PROPAGATION AND 
INSTABILITIES: II. OBLIQUE PROPAGATION 

by 
t 

J. A. Tataronis and F. W. Crawford 
Institute for Plasma Research 
Stanford University 
Stanford , California 

ABSTRACT 


Part I of this paper has discussed longitudinal wave propagation, 
perpendicular to the static magnetic field, in a warm homogeneous 
magnetoplasma. Part II extends the work to oblique propagation. The 
ring and Maxwellian transverse velocity distributions, which are 
absolutely unstable and stable, respectively, for perpendicular prop- 
agation, are shown here to be both absolutely unstable for oblique 
propagation. Their instabilities are radically affected by the intro^ 
duction of axial velocity spread. A detailed study is made for the 
cases of a Maxwellian transverse distribution with resonance and 
resonance-squared axial velocity distributions. It is shown that 
absolute instabilities occur at high energy anisotropies. These become 
convectively unstable as the anisotropy is reduced, and are quenched 
as isotropy is approached. The isotropic Maxwellian is treated in 
detail, and shows both cyclotron and Landau collisionless damping. 
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1. IOT’RODUCT ION 


In Part I (Tataronis and Crawford 1969), we considered the prop- 
erties of perpendicularly propagating CHW for the ring, spherical shell, 
and Maxwellian electron velocity distributions, and for a mixture of 
ring and Maxwellian. We continue by considering oblique propagation. 

In Section 3, a thorough study of the ring distribution is made. This 
is followed by consideration in Section 4 of a Maxwellian transverse 
distribution. Though highly instrictive, these distributions are not 
realistic in the sense that the electrons in practical magnetoplasmas 
will always have some motions parallel to the static magnetic field. 

For this reason, the remainder of the paper is devoted to assessing 
the influence of such motions on the strength and nature of CHW 
instabilities. Section 5 retains the Maxwellian transverse distribution, 
and investigates the effects introduced by resonance and resonance- 
squared axial velocity distributions. This requires some refinements 
to the stability criterion used in Part I. These are presented in 
Section 2. Finally, Section 6 treats an isotropic Maxwellian distribution 
As before, we emphasize that we are not presenting a review of 
previous work, but are trying to present systematically in a single 
source the salient features emerging from our numerical studies of the 
very complicated CHW dispersion relations appropriate to a variety of 
basic electron velocity distributions. The pioneer work in the field, 
and many more recent results, should be accessible via the review by 
Crawfo rd (1967) . 



2 . THEORY 


The general CHW dispersion relation may be written as (see Part I ) 



For velocity distributions of the form f^(v ,v ) = f (v )S(v ). 

O' 1> H 7 i' l' H 7 

( 1 ) reduces to 



( 2 ) 


Two such cases, i.e. without axial velocity spread, will be studied in 
Sections J and k. It is clear from ( 2 ) that the roots <r>(.k real) must 
be either real, or occur in complex conjugate pairs. This is similar to 
the behavior of the cases treated in Part I. Mode coupling, and absolute 
instability, will again be shown to occur. 

When there is axial velocity spread, the situation is more compli- 
cated, and it becomes important to note that K(^,k) given by (l) is not 
a unique function, but has two branches for k real. These are readily 

’ 11 

separated by introduction of the identity 
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where = {&> - nn> c )/k^ . To ensure convergence of the integral, the plus 
sign is chosen if w_ < 0, and the negative if w^ > 0. Substituting 
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( 3 ) in (l), and assuming cd^ < o, leads to 
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Here, H ( C,k ,k ) is the Fourier transform. of H (v ,k ,k ) , 
’ n k l’ || rr || J l’ || 


( 4 ) 
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H (£,k k ) = / dv exp( iv £) H ( v ,k k ) . (5) 

n l ’ II 7 II II n II 1 II 


. -p 

The two branches K ^ K , are connected by the relation 


K + (-d>* k ,k ) = K~(o> k -k )* . 

; X’ II 1 || 


( 6 ) 


As a. result of the decomposition of K(oo^k ,k ) in ( k) , the insta- 

1 It 

bility criterion stated in Part I, Section 2, must be modified. The appro 

priate changes have been stated by Derfler ( 1967 ) . For discussion , assume 

that k^ is a fixed, real number and that the propagation is parallel to 

the magnetic field, Derfler* s work then implies that absolute instability 

is present if zeros of K + (<i),k ,k ) collide across the positive k axis 

’ -i-' ii it 

or if zeros of K (0),k ,k ) collide across the negative k axis, as 

' ii ii 

varies along some contour in the lower half complex plane. We shall make 
use of this modified instability criterion in later sections. 
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3. RING DISTRIBUTION 


This distribution may be written as (see Part l) 


f o ( vO 


5 5(v •-* v ) 8 (v ) 

2ltv 0i 1 0x II 


(7) 


Substitution of (j) in ( 2 ) yields the dispersion relation 
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Remembering that k = k^ + k } ( 8 ) may be rearranged to give 
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where 0 is the angle between k and the static magnetic field. 

Dispersion Characteristics 

Figure 1 shows numerical solutions of ( 19 ) for several values of 0, 

2 2 o 

and CD /CD = 1. The curves for 0 / Q0 indicate that there are now two 
pc 

branches undulating about each cyclotron harmonic line. At ncD , there 

c 

are resonances (all n), and cutoffs ( n / 0) . Two additional cutoffs can 
be located by putting ji = 0 in ( 9 ), 

CD 2 = i ! CD 2 + CD 2 ± [(cd 2 -cd 2 ) 2 + ijCD 2 CD 2 sin 2 0]^l . (lO) 

2 ( p c LV p c pc ) ' ' 

2 2 “ 

These tend to (£ (& as 6 tends to zero, and (co +ar ) 2 zero, as 9 

p* c 7 7 pc 7 7 

tends to 90 ° . 

It is easy to see from ( 9 ) that the undulating branches cut the har- 
monic lines at the zeros of J (u ) . As Q tends to 90° 

n jl 


9 


the curves 




Fig. 1, DISPERSION CHARACTERISTICS OF OBLIQUELY PROPAGATING 
CYCLOTRON HARMONIC WAVES: RING DISTRIBUTION. 
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pinch towards the points J^(p. ) = 0, so that in the limit one undulat- 
ing branch intersects each harmonic line at J (u ) = 0 and J^u. ) = 0. 

n'*x n' r i 

We note from the curves for 9 = k-5° that a downgoing loop from 0)/to = 1 

c 

has coupled with an upgoing loop from =0 to produce an absolute 

c 

instability. An additional instability is exhibited for 9 - 15° , where 
for convenience (dashed) has been plotted using tD/^ c = 1 as a 

base line. We conclude from these results, and table 1 of Part 1, that 
instability may occur in oblique propagation well below the threshold for 
unstable perpendicular propagation. 

A thorough study of ( 9 ) could be made via. diagrams such as figure 1. 
Since we are solving for Ch(k real), and the instabilities are absolute, 
we have preferred the alternative of producing the families of curves 
shown in figure 2 for fixed k^ and varying .k > These have the advan- 
tages of being easily related to the work of Part I, and are directly 

relevant to parallel propagation in a transversely bounded plasma, where 
the finite dimensions could be accounted for roughly by taking a fixed 
value of k^. 

To find the cutoffs for figure 2, we put k = 0 in ( 9 ). One set 
of cutoffs are given by the zeros of the numerator, which is just the 
CHW dispersion relation for perpendicular propagation, so they may be 
read off figure 2 of Part I. Our choice of u = k v /co = 1 , 3 4.5, 

will now be appreciated. The first value is in a region where no insta- 

bility can occur for perpendicular propagation; the second is where only 
the zero-frequency instability can occur ^ and the third is where mode 
coupling can occur for 1 ^ 03 ^/cd^ ^ 2. A second set of cutoffs can be 
obtained by approximating ( 9 ) by 



2n J / (p. ) 
n X 




(11) 


Thus cd/cd -» n as k -♦ 0, so there is a cutoff at each cyclotron 
' c || * 

harmonic. 

The resonances (k = oo) occur when the denominator of ( 9 ) is zero. 

II 2 2 

This function is plotted in figure 3 a - For 03 r>n ~ , two branches 

P C 
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Fig. 2. CONTINUED. 
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Fig. 2. CONTINUED. 
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(a) Ring distribution (b) Maxwellian transverse 

distribution 


Fig. 3. RESONANCE FREQUENCIES (k |( = «>) OF OBLIQUELY 
PROPAGATING CYCLOTRON HARMONIC WAVES. 
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undulate about each harmonic line^ and intersect them at the zeros 

(a ) of J (p. ). There is no coupling between upgoing and downgoing 
nm n op 

loops, so the resonance frequencies are real. For = 1, coupling 

occurs. For all cases we have considered, coupling occurs in the ranges 

and a, , .. N < p < a , , x , f or <o 2 /<n 2 suf f i- 

V n+l),m **jl o n,(m+l) * p c 


a < u < a 

nm X 1 


( n-f 1 ) 7 m 

ciently large* Some critical values of 


CD' 


V62 


are given in table 1 for 


the first three frequency bands. They are clearly much lower than those 
for purely perpendicular propagation (see table 1 of Part i). 


Frequency Band 


r 0 < 03/0i < 

c 

1 

1 < <x>/<x> < 2 

c 

2 < CD/(b < 3 

C ^ 


Range of 

p p 

03 /CD 
P c 

Range of 

p p 

03 VCD 
P C 

Range of ll 0 > 2 /cd 2 

^X pc 

0 - 2.40 

0.34 

0 - 3.83 

0.53 

0 - 5.14 

0.66 

2.40 - 5.83 

1.07 

3.83 - 5.14 

1.70 

5.14 - 6.38 

2.29 

3.83 ~ 5.52 

1.38 

5.14 - 7.02 

1.62 

6.38 - 8.42 

1.82 


TABLE 1 . Instability threshold values of cd^/cd^ 


We now consider the families of dispersion characteristics displayed 

2 2 

in figure 2. For cd /cd small, each passband has two propagation 

pc 7 

branches associated with it,, the frequency spread about no3 c decreasing 
rapidly with increasing n. For jjl^ = 1, instability first occurs in 
the lowest passband , for a value of k v q^° c corresponding to propaga- 
tion at 0 ~ 45° . Successive passbands become unstable roughly for 
2 2 2 

p + *c > (^c } 9 S rowt3:1 ra * tes are small except in the first 

few passbands. For j i ~ 3, instability occurs first in the second and 

-1 pp 

third passbands, and for 03 /03 = 3 has spread to the first and fourth. 

22 p c 

Above .(13 /cd as 17,02, we expect additional effects since the zero- 
p c 7 

frequency instability has set in for purely perpendicular propagation. 

2 2 

These appear in the lowest passband for cd^/cd^ 5=5 20. At small 

there are two imaginary roots. As k increases ^ there is a transition 

to complex roots with a very rapid growth rate. Further increase in 

2 2 

cd /CD reduces the real parts to zero, leaving purely imaginary solutions 
pc 

for all k . For the final case, u. = 4.5, where unstable perpendicular 

11 

propagation first occurs in the band 1 ^ <d /to ^ 2, the characteristics 
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in all other bands are similar to those for \± ^ == 1, though with stronger 
growth rates. In the second pa.ssband, however, there is absolute insta- 
bility for all k once the perpendicular instability threshold has been 

II 2 2 

passed. The results for /cd ^8 suggest that the instability in this 

p ^ 

band always has the highest growth rate. 

A common feature of the curves in figure 2 is that the branches 
0)(k^) increase or decrease monotonically towards coupling points. It 
follows that the threshold for coupling will have k = oo, and that 
curves such as those of figure 3 determine immediately whether a given 
band will be unstable or not. The instability thresholds suggested in 
this way are substantially lower than for perpendicular propagation. In 
anticipation of Section h y it should be pointed out that, in practice, 
there is always a spread of electron velocities parallel to the magnetic 
field. This can strongly modify the threshold conditions. An additional 
caution should be stated against interpreting too literally the limit 
k -* oo since the theory will ultimately break down at wavelengths com- 

11 

parable with an electronic Debye length. 
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k. MAXWELLIAN TRANSVERSE DISTRIBUTION: NO AXIAL VELOCITY 


In Part I, it was established that introduction of a transverse 
velocity spread led to less violent perpendicularly propagating insta- 
bility than that associated with a ring distribution, and could give 
complete stability if hf Q /hv^ < 0 for all v^ >0. The Maxwellian 
falls into this class, and we now wish to determine whether stability 
persists when oblique propagation is considered. 

For no axial velocity spread, the distribution may be written as 


f (v ,v ) = 
1 l’ II 


exp - 


2 tcv' 


ti 


2 v‘ 


2 f ( V 

tj. , 


( 12 ) 


Substitution in (2), and some manipulation, yields the dispersion rela- 
tion 


2 oo 
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2 2 
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nciD (co - ntD Y 
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where \ = (k v /<& ) . Equation ( 13 ) may be rearranged analogously to 
X t) x c 

(9) as 


03 zL 


cot^0 = = 


0) n=-oo 
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03 - n03 


k 


a ) 2 00 
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X '^2 2 a!P< - X)l » (X) 


CD 
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Dispersion Characteristics 

As for the ring distribution, one set of cutoffs are found from (.1.4) 

at no) , and a second set are determined by the numerator. Its zeros 
c’ 

are solutions of the CHW dispersion relation for perpendicular propaga- 
tion, and can be read off figure 5 of Part I. Unlike the ring distribu- 
tion, these cutoff frequencies are always real. Resonances occur at 
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zeros of the numerator of (l4). Their behavior is shown in figure 3 b, 
contrasted with that of the ring distribution. There are two modes in 
each passband, but instead of undulating about the harmonic lines, their 

separation increases to a maximum then decreases to zero as k -» °° . 

2 2 

Even for cd/cd = 1 , coupling is possible, and figure ya shows absolute 
instability in the first two passbands. A consequence of the nonundula^ 
toiy nature of the propagation branches is that there can be instability 
in only a single range of k v /cd , for a given passband. For distri- 
butions such as the ring, there may be many ranges. 

Figure 4 shows typical dispersion characteristics computed from 
(l4) . They may be compared with the relevant curves of figure 2b for 
the ring distribution, and are found to be very similar. Although an 

exhaustive numerical study has not been made, we can be sure that there 

2 2 

are no long wavelength (k — 0 ) instabilities for any value of CD /CD. . 

II pc 

Our next concern must be to see to what extent the short wavelength 
instabilities are quenched by the introduction of axial velocity spread. 
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Fig. 4. DISPERSION CHARACTERISTICS OF OBLIQUELY PROPAGATING 
CYCLOTRON HARMONIC WAVES: MAXWELLIAN TRANSVERSE DISTRIBUTION. 

NO AXIAL VELOCITY. 
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5. MAXWELLIAN TRANSVERSE DISTRIBUTION: RESONANCE AXIAL DISTRIBUTION 

To derive the CHW dispersion relation in the presence of a velocity 
spread parallel to the static magnetic field, (l) is required, and for 
axial distributions such as the Maxwellian leads to an expression trans- 
cendental in k^. The numerical stability analysis is then greatly com- 
plicated by the presence of an infinite series of Landau poles. These 
mathematical difficulties can be avoided, with little loss to the phy- 
sics, by considering plausible analytic approximations to the axial velo- 
city distribution of interest, constructed so as to yield simpler disper- 
sion relations algebraic in k , i.e. with a small number of Landau 

II 

poles. For this reason, we shall now investigate resonance-type distri- 
butions of the forms 


f n( v ,> v J = 
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Substitution of the first of these in (l) yields 


or 
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Substitution of the second ("resonance-squared") distribution in (l) 
leads to 


or 
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Stability Analysis 


Comparison of (l6) with (13) shows that the dispersion relation for 
the simple resonance axial distribution can be obtained directly from 
that with no axial velocity by replacing cu with a>-ik v . Since the 

, Jl H , 

initial step in the stability analysis is to solve (l6) for oo(k real.), 
we see that the transformation changes o^, but not a> . Figure 5 
illustrates this for the parameters of figure 4 . Both of the imaginary 
parts of the complex solutions are plotted. The real parts are still 
given by figure 4 . We note that the introduction of the complex Doppler 
shift due to axial velocities has the effect of stabilizing the short 
wavelengths (k -*<»). In the first passband, there is now only a narrow 
range of k^v^/o^ for which cd_^ < 0 . In the second passband, the range 
is even more restricted. The third passband has been stabilized, and all 
higher passbands (3 < show the same attenuation, = k^v^. 

In previous sections, and Part I of the paper, complex conjugate 
solutions were obtained for o>(k real), and the instabilities were inter- 

cw * 

pretable as absolute. With an axial velocity spread, this is no longer 
the case ; and more detailed stability analysis is required. This could 
be carried out for the simple resonance distribution, but we shall not 
do so, since it has the undesirable feature of an infinite thermal velo- 
city. We have employed it because of the direct demonstration it pro- 
vides^ and the insight it gives, into the stabilizing effect of axial 
velocities. For the stability analysis, we shall consider the dispersion 
relation of (17), for the resonance-squared distribution. This has a 

thermal velocity of v . 

t ll 

Figure 6 shows typical dispersion characteristics computed from (17) 
for the parameters of figure k. As in our previous examples, there are 
two propagating branches in each passband. For clarity, these have been 
separated into two sets. We note that there is one unstable branch in 



and a root with dipping down towards instability in the second 

1 c 

passband. 


The latter is studied in more detail in the conformal mappings of 
figure 7 , which show contours A-D in the complex a)- plane mapped into the 
complex k -plane via (17). In accordance with the stability criterion 
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k 1l , *l / “e 


Fig. 5. DISPERSION CHARACTERISTICS OF OBLIQUELY PROPAGATING 
CYCLOTRON HARMONIC WAVES: MAXWELLIAN TRANSVERSE DISTRIBUTION. 

RESONANCE AXIAL DISTRIBUTION. 
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Fig. 6, DISPERSION CHARACTERISTICS OF OBLIQUELY PROPAGATING 
CYCLOTRON HARMONIC WAVES: MAXWELLIAN TRANSVERSE DISTRIBUTION. 

RESONANCE-SQUARED AXIAL DISTRIBUTION. 
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<*> Vv = °- 42 (b) V v u - °- 43 (c> v t«/ v u = °- 50 


Fig. 7. STABILITY ANALYSIS OF OBLIQUELY PROPAGATING CYCLOTRON 
HARMONIC WAVES: MAXWELLIAN TRANSVERSE DISTRIBUTION. RESONANCE- 

SQUARED AXIAL DISTRIBUTION. 
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quoted in Section 2, we look for saddlepoints (BK + /Bk = O) correspond- 

" 4 ~ / \ ^ 

ing to the merging of two zeros of K (<*> k) across the positive k axis, 

II 

and their corresponding branchpoints in the cu-plane. For v /v = Q.k2 

t \\ tA 

merging occurs with 03. < 0 at the branchpoint, so there is absolute 

1 7 

instability. An increase to v, /v A = 0.^3 makes co. >0 at the 

1 1| tx v 1 

branchpoint. Since there are still solutions with < 0 for 

k . v. M = 0, the situation is now convectively unstable. For a fur- 
||X t|| c > 

ther increase, to v, /v, = 0 . 50 , we see that W. >0 for k . v /to =0 

' A II tJ - 1 II 1 All c 

This implies that the convective instability has been quenched and the 


root is stable. There would still remain of course, the unstable root 


shown in the first passband in figure 6. This could be removed by a fur- 
ther increase in v, /v, . 

A II AX 

To complete the stability analysis, it would be necessary to fix 

v , /v, and repeat the procedure for all real k . We will not do this, 

1 1| tx i 9 

but will simply state the following speculative general conclusions, 
which are supported by other computations we have carried out on this 
and other magnetoplasma instabilities: Oblique propagation always leads 

to absolute instability if there is zero axial velocity spread, i.e, 
infinite anisotropy. As the anisotropy is reduced, the instabilities 
become convective. They are finally quenched as isotropy is approached. 
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6 . ISOTROPIC MAXWELLIAN DISTRIBUTION 


To complete the comparison with the basic cases treated in Part I, 
we now consider a Maxwellian transverse velocity distribution, and a 
Maxwellian axial distribution. Most generally, this may be written as 




(at) 


3/2 2 

v . v, 
t u ti 


expl- 


2v 


tJL 



( 18 ) 


where v and v are the parallel and perpendicular thermal veloci- 

t|j tj. 

ties, respectively. Substitution in (1) leads to the dispersion relation 




K(co,k) = 1 + 


.2 2 

kV t, 


00 


1 + ^ exp(-x)l n (x) 


to - (l-T)n£D /cn - ncD 


Z| 


n=-oo 


22k V 


= 0 


(19) 


2 2 

where T = v^ /v, , and Z(z), which arises from the v integration, 

tn tx> , II 

is defined as (Fried and Conte 1961) 


00 . / cr=0 , z ± < 0\ 

Z(z) = -^ f dt - Ict rt2 exp(-z 2 ) . =1 , = 0> . (20) 

**-00 " Z =2, >0) 

Here, the principal part of the integral is taken if z lies on the 
real t axis. 

For anisotropic plasmas with v < v , instabilities similar to 

L || TJL 

those discussed in Section 5 are to be expected. They have attracted 

considerable interest , and the reader is referred to the literature 

(see Crawford 1968 for references). Here we shall consider the isotropic 

case (v^ = = v . ) . which is completely stable (Bernstein 1958). 

t t|j tx 7 

Dispersion Characteristics 

Figure 8 shows numerical solutions o>(lc real) to (19), and con- 
firms the absence of instability. As in all other cases in Part II , 
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0.25 0.50 0 0.25 0.50 0 0.25 0.50 0 0.25 0.50 0 0.25 0.50 

k ll v t /w c 


Fig. 8. DISPERSION CURVES FOR OBLIQUELY PROPAGATING CYCLOTRON 
HARMONIC WAVES: ISOTROPIC MAXWELLIAN DISTRIBUTION. 



there are two branches in each passband. These have been separated, for 
clarity. Although there is no collisionless damping for purely perpendi- 
cular propagation, we note that the attenuation increases very rapidly 

with k v,/Cb particularly for the modes with cutoffs at m) . Bear- 
11 t c ; 17 c 

ing in mind that 0) corresponds to an attenuation of over 50 dB/ 

cyclotron period (= 2rt/<^ c ), it will be appreciated from figure 8 that 
propagation will be difficult to observe experimentally for angles only 
a few degrees off exact perpendicularity. 

Since there are no absolute instabilities present, we are entitled 
to solve ( 19 ) for k(<+> real). The character of the Z-function of (20) 
makes this by no means trivial, however and for complex k we shall pre- 
sent only the set of curves shown in figure 9- Particular branches have 
been chosen so as to illustrate the passage from no collisionless damping 
in the limit 0 = 90°, to Landau damping of plasma oscillations when 
0=0. A distinction may be made in this respect between cyclotron and 
Landau-type interactions: An electron spiralling about the magnetic 

field sees the wave at the Doppler-shifted frequency cu - k v^. If a 

spread in v exists, it is possible to find electrons for which 

I! 

cd ! na) c ; where n is an integer. There is then an energy exchange 
between the electrons and the CHW via its electric field. Damping 
results if the electrons gain energy, while instability occurs if the 
electrons lose energy. Cyclotron damping (or instability) arises from 
those electrons which see the CHW at no> (n 4 0). If the electrons 
see a static field, corresponding to n = 0, Landau damping (or insta- 
bility) results. In the limiting case 0=0, the electric field of 
the wave is parallel to B_ . so the damping is purely of the Landau 
type. 
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7, DISCUSSION 


The computations of Sections 3-6 demonstrate that oblique prop- 
agation has features which disappear in the perpendicular propagation 
limit. Most important is the presence of a second branch in each 
cyclotron harmonic band. These contribute materially to mode coupling 
and the occurrence of instabilities. For the ring distribution, the 

results of Section 3 indicate much lower instability thresholds 
2 2 

(a) /cu ), while Section 4 demonstrates that instability can occur in 
p c 

oblique propagation for the Maxwellian transverse distribution; a 
distribution which was shown in Part I to be stable for perpendicular 
propagation. 

Addition of an axial velocity distribution produces successively 
three new effects. First, it stabilizes short wavelength instabilities. 
Second, it can convert absolute to convective instability, and third 
it can quench instabilities entirely if isotropy is approached. Indeed, 
the dispersion relation for an .isotropic Maxwellian generally shows 
very strong damping, except within a few degrees of exactly perpendicular 
propagation. 

In Part I, a mixture of ring and Maxwellian distributions was 
examined in detail . It would be very time-consuming and expensive to 
extend this work to oblique propagation, but we may speculate on the 
likely behavior of this and other mixtures of energetic electrons in a 
stable background. It seems likely that the high damping in a Maxwel- 
lian plasma could only be offset by a relatively dense energetic group, 
and that marginal stability would occur for nearly perpendicular prop- 
agation. Further work is required, and may point the way to a suitable 
experiment for checking- the theory. For example, the density of the 
energetic group could be raised while oblique GHW propagation was 
measured in a Maxwellian background plasma. The reduction of damping 
and the onsets of convective and absolute instability should all be 
measurable for comparison with numerical predictions. 

This work was supported by the U.S . Atomic Energy Commission, and 
the National Aeronautics and Space Administration. The authors are 
indebted to Dr. H. Derfler for many helpful discussions. 
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